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ANALYSIS  OF  STEERABILITY  OF  TRACKED  VEHICLES; 
THEORETICAL  PREDICTIONS  VERSUS  FIELD  MEASUREMENTS 


PART  I:  INTRODUCTION 

1.  Development  of  high-mobility/agility  tracked  combat  vehicles  has 
received  considerable  attention  recently  because  of  the  possibilities  they 
offer  for  increased  battlefield  survivability  through  the  avoidance,  by 
high-speed  and  violent  maneuver,  of  hits  by  high-velocity  projectiles  and 
missiles.  In  order  to  design  and  develop  such  vehicles  rationally,  it  is 
necessary  to  have  a  quantitative  understanding  of  the  interrelationship 
between  the  terrain  factors  (such  as  soil  type,  soil  shear  strength  and 
compressibility,  etc.)  and  the  vehicle's  characteristics  (weight,  track 
length  and  width,  location  of  center  of  gravity,  velocity,  etc, )  during 
steering.  To  study  such  an  interrelationship,  it  is  necessary  to  construct 
idealized  mathematical  models  of  the  terrain-vehicle  interaction.  The 
accuracy  and  range  of  application  of  such  models  must,  of  course,  be  deter¬ 
mined  from  actual  mobility  experiments  and  obviously  must  depend  on  the 
degree  of  relevance  of  the  idealized  model  as  an  approximation  to  the  real 
behavior. 

2.  The  basic  concepts  of  the  theory  of  terrain-vehicle  interaction 
were  developed  by  Bekker  during  the  1950’ s  (Reference  1).  By  assuming 
various  load  distributions  along  the  tracks,  Bekker  was  able  to  develop 
several  mathematical  expressions  relating  the  characteristics  of  the  vehicle 
and  the  tractive  effort  of  the  terrain  during  steering.  By  considering  the 
lateral  and  longitudinal  coefficients  of  friction  between  the  track  and  the 
ground,  Hayashi  (Reference  2)  developed  simple  equations  for  practical 
analysis  of  steering  of  tracked  vehicles.  Hayashi' s  work,  however,  did  not 
include  the  effect  of  the  centrifugal  forces  on  steering  performance  of  the 
vehicle.  Kitano  and  Jyorzaki  (Reference  3)  developed  a  more  comprehensive 
model  for  uniform  turning  motion  including  the  effects  of  centrifugal  forces. 
This  model,  however,  is  based  on  the  assumption  that  ground  pressure  is 
concentrated  under  each  road  wheel  and  the  terrain-track  interaction  is 
simulated  by  Coulomb-type  friction.  The  model  given  by  Kitano  and  Jyorzaki 
was  extended  by  Kitano  and  Kuma  (Reference  4)  to  Include  nonuniform  (tran¬ 
sient)  motion,  but  the  basic  elements  of  the  terrain-track  interaction  part 


of  the  model  were  retained,  iialedi  and  Rohani  (Reference  5)  developed  a 
model  for  uniform  turning  motion  pa  allel  to  the  development  reported. in 
Reference  3  insofar  as  the  kinem.jtics  of  the  vehicle  are  concerned.  In 
contrast  to  Beference  3,  however,  this  model  is  based  on  a  more  comprehensive 
soil  model.  1*0  the  present  q^peir,  thn  terrain-vehicle  model  reported  in 
Reference  5  is  extended  to  include  n  jiiuniform  (transient)*  motion.  In  addi¬ 
tion,,  the  soil  model  is  modified  to  Include  a  nonlinear  failure  envelope 
describing  the  shearing  strength  of  the  te*  rain  material. 

3.  To  demonstrate  the  application  of  the  model,  the  steering  perform¬ 
ance  of  an  armored  personnel  carrier  has  b.--'.en  predicted  and  correlated  with 
full-scale  test  results. 
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PART  II:  SOIL  MODEL 


4.  One  of  the  most  important  properties  of  soil  affecting  traffic- 
ability  is  the  in  situ  shear  strength  of  the  soil.  The  shear  strength  of 
earth  materials  varies  greatly  for  different  types  of  soil  and  is  dependent 
on  the  confining  pressure  and  time  rate  of  loading  (shearing).  This  depen¬ 
dence,  however,  is  not  the  same  for  all  soils  and  varies  with  respect  to  two 
fundamental  strength  properties  of  soil;  the  cohesive  and  the  frictional 
properties.  It  has  been  found  experimentally  that  the  shear  strength  of 
purely  cohesive  soils  (soils  without  frictional  strength)  is  independent  of 
the  confining  stress  and  is  strongly  affected  by  the  time  rate  of  shearing. 

On  the  other  hand,  in  the  case  of  purely  frictional  soil  (soils  without 
cohesive  strength),  the  shear  strength  is  found  to  be  independent  of  time 
rate  of  loading  and  is  strongly  dependent  on  the  confining  pressure.  In 
nature,  most  soils  exhibit  shearing  resistance  due  to  both  the  frictional  and 
cohesive  components.  The  cohesive  and  frictional  components  of  strength  are 
usually  added  together  in  order  to  obtain  the  total  shear  strength  of  the 
material;  i.e.. 


=  A  -  M  exp  (-No) 


(1) 


where  is  the  maximum  shearing  strength  of  the  material,  C  =  A  -  M  is 

the  cohesive  strength  of  the  material  corresponding  to  static  loading  (very 
slow  rate  of  deformation),  a  is  normal  stress,  and  N  is  a  material 
constant.  Equation  1  is  shown  graphically  in  Figure  1. 


Effect  of  Rate  of  Deformation 

5.  As  was  pointed  out  previously,  the  cohesive  strength  of  the  mate¬ 
rial  is  dependent  on  the  time  rate  of  loading  (shearing);  i.e.,  the  cohesive 
component  of  strength  increases  with  increasing  rate  of  loading.  For  the 
range  of  loading  rates  associated  with  the  motion  of  tr^icked  vehicles,  the 
contribution  to  cohesive  strength  due  to  dynamic  loading  can  be  expressed  as 
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Figure  2.  Proposed  soil 
stress/deformation  relation 
during  shearing  process. 
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C^£l  -  ei:p(-AA)2  ,  where  and  A  are  material  constants,  and  A  is 

time  rate  of  shearing  deformation.  In  view  of  the  above  expression  and 
Equation  1,  the  dynamic  failure  criterion  takes  the  following  form: 

« 

Tj^  =  A  +  C^[l  -  exp(-AA)]  -  M  exp  (-No)  (2) 

Shear  Stress-Shear  Deformaton  Relation 


6.  Prior  to  failure,  the  shear  stress-shear  deformation  character¬ 
istics  of  a  variety  of  soils  can  be  expressed  by  the  following  mathematical 
expression  (Reference  6) : 


(3) 


The  behavior  of  Equation  3  is  shown  graphically  in  Figure  2.  In  Figure  2, 

T  denotes  shearing  stress,  A  is  shearing  deformation,  and  G  is  the 
initial  shear  stiffness  coefficient.  In  view  of  Equation  2,  the  shear 
stress-shear  deformation  relation  for  soil  .(Equation  3)  becomes 

G[A  +  C ,  -  C ,  exp(-AA)  -  M  exp(-No)jA 

T  - - - - - - ^ -  (4) 

G|a|  +  a  +  -  Cj  exp(-AA)  -  M  exp (-Wo) 


For  purely  cohesive  soils,  N  =  0  and  t  is  only  a  function  of  A  and 
A  .  For  granular  material,  M  =  A  and  is  zero,  and  t  is  a  function 

of  A  and  a  .  For  mixed  soils  exhibiting  shearing  resistance  due  to  both 
frictional  and  cohesive  components,  t  is  dependent  on  A  ,  A  ,  and  a  . 

An  appropriate  test  for  determining  the  numerical  values  of  the  six  material 
constants  in  Equation  4  is  an  in  situ  direct  shear  test.  A  field  direct 
shear  device  has  been  developed  at  WES  for  this  purpose.  A  description  of 
this  device  and  the  method  of  analysis  of  the  data  obtained  from  the  direct 
shear  test  are  documented  in  Appendix  C  of  Reference  7.  In  the  following 
section,  the  equations  of  motions  for  a  track-laying  vehicle  during  steering 
are  developed  using  the  proposed  soil  model  (Equation  4)  in  conjunction  with 
track  slippage,  centrifi^gal  forces,  and  vehicle  characteristics. 
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PART  III:  DERIVATION  OF  TERRAIN-VEHICLE  MODEL 

\ 

Boundary  Conditions 

7.  The  geometry  of  the  vehicle  and  the  boundary  conditions  of  the 
proposed  nodel  are  shown  schematically  in  Figure  3.  The  XYZ  coordinates  are 
the  local  coordinate  system  of  which  X  is  always  the  longitudinal  axis  of 
the  vehicle  and  Y  is  a  transverse  axis  parallel  to  the  ground.  These  axes 
Intersect  at  the  center  of  geometry  of  the  vehicle  0.  The  Z  axis  is  a 
vertical  axis  passing  through  the  origin  0.  The  center  of  gravity  of  the 
vehicle  (CG)  lies  on  the  X  axis  and  is  displaced  by  a  distance  C^^  from  the 
origin.  The  numerical  value  of  C^^  is  assumed  to  be  positive  if  CG  is 
displaced  forward  from  the  center  of  geometry  of  the  vehicle.  The  XY 
coordinates  of  the  instantaneous  center  of  rotation  (ICR)  are  P  +  C^  and 

A 

R  ,  respectively,  where  P  is  the  offset-  The  center  of  rotation  and  the 
radius  of  the  trajectory  of  the  CG  are,  respectively,  CR  and  R^.  The  height 
of  the  center  of  gravity  measured  from  ground  surface  is  denoted  by  H.  The 
lengths  of  the  track-ground  contact,  the  track  width,  and  the  tread  of  the 
tracks  are  L,  D,  and  B,  respectively.  As  shown  in  Figure  3,  the  components 
of  the  inertial  forces  F„  in  X  and  Y  directions  are,  respectively,  F„,,  and 
F^^.  The  weight  of  the  vehicle  is  W. 

'  Stress  Distribution  Along  the  Tracks 

8.  Two  types  of  stress  (i.e.,  normal  and  shear  stresses)  exi«st  along 
the  track.  As  indicated  in  Figure  3,  the  normal  stresses  under  the  outer 
and  inner  tracks  are  denoted  by  R^(X)  and  R2(X),  respectively.  The  compo¬ 
nents  of  the  shear  stress  in  X  and  Y  directions  are,  respectively,  Tj^(X) 
and  Qj^(X)  for  the  outer  track,  and  T2(X)  and  Q2(X)  for  the  inner  track. 

These  stresses  are  dependent  on  the  terrain  type,  vehicle  configuration,  and 
speed  and  turning  radius  of  the  vehicle. 

9.  The  magnitude  of  normal- stresses  R^(X)  and  R2(X)  can  be  determined 
in  terms  of  the  components  of  the  inertial  force,  the  track  tensions,  and 
the  characteristics  of  the  vehicle  by  considering  the  balance  of  vertical 
stresses  and  their  moments  in  Figure  3.*  Thus: 

*  For  sake  of  brevity,  the  effect  of  track  tension  is  not  included  in  this 
paper  and  was  "turned  off"  in  the  computer  program  AGIL.  The  reader  is 
referred  to  Reference  7  for  a  complete  analysis  of  track  tension  and  its 
effect  on  steering  performance 'of  tracked  vehicles. 
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where  h  ,"  H/L  ,  b-»  B/L  ,  d-  *  D/L  ,  =  C^^/L  j  x  =  X/L  ,  y  =  Y/L  ,  and 

z  *  Z/L  . 

10.  The'-!components-of;.ithe  shear  stress  in-  the  X‘ atid  ’Y  directions 
along.‘both  the -outer  and  inner- tracks  can-'be  obtained  by  combining;  Equations 
4,  5,  and  6.*  Thus  (it  is  . noted  that  and,  R2i  replace  the  normalistress  a 
.in  Equation  44)  :=  : :  .  , 


Ij(x)  -  ^ 


da+dc^-dcj  exp(-X6£)-m  exp[-nrj|^(K)] 


cos  Yx 


L  (p|6j^|d+da+dCj*?dCj  exp(-X6^)-m  exp[-nr^(x)] ' 


„  da+dc.-dc,  exp(-X6. ).-m  exp[~nr .  (x)] 

Q..(x)  -  \  p6  - r- - S - S - L_ - i - sin  Y.  (8] 

L  jj I  6 I d+da+dc^-dcj  exp(-X6j^)-m  exp[rnr^(x)3 

where  i  «  1,2  ;  rj.(x)  »  dL^R^(x)/W  ;  6^  ■  A^/L  ;  6^,  ■  ;  y  =  GL^/W  ; 

X  *  AL  5  a  =  AL^/W  ;  m  -  ML^/W  ;  n  -  NW/L^  ;  and  c^  -  C^L^/W  .  The 
variables  Yj^  and  v.^  »  Equations  7  and  8,  are  the  slip  angles  and  can 
be  written  as 


X  -  P  -  C 


=  tan 


’X  -I’^-P-.'^X 


=  tan 


-1  *  -  ^  -  ‘=x  -1  *  -  P  -  'x 

I'a  ■  - -  - - 


where  =  Cj^/L  ,  ~  ^2^^  *  P  “  parameter  is  the 

distance  between  the  Instantaneous  center  of  rotation  of  the  outer  track  and 
its  axis  of  symmetry,  and  C2  is  the  distance  between  the  instantaneous 
center  of  rotation  of  the  inner  track  and  its  axis  of  symmetry. 


*  To  account  for  the  effect  of  the  size  of  the  shear  box  on  the  shear  stiff¬ 
ness  G  ,  the  measured  value  of  G  is  normalized  in  AGIL  by  multiplying  it 
by  4/L  (the  length  of  the  shear  box  =*  4  in.). 
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11.  In  order  to  use  Equations  7  through  9,  the  track  slip  velocities 

•  « 

and  displacements  (i.e.,  ,  A^^  ,  j  and  and  the  inertial 

forces,  and  have  to  be  determined. 

Kinematics  of  the  Vehicle 


12.  A  tracked  vehicle  in  transient  motion  is  shown  schematically  in 

Figure  h.  The  XYZ  coordinates  are  the  local  coordinate  systems  that  are 
fixed  with  respect  to  the  moving  vehicle  (also  see  Figure  3).  The  origin  0 
of  this  coordinate  system  stays,  for  all  time,  at  a  distance  from  the 
center  of  gravity  of  the  vehicle.  The  ^coordinate  system  is  fixed  on 
level  ground,  and  its  origin  coincides  with  the  center  of  gravity  at  time 
zero.  The  vehicle  can  maneuver  on  the  plane  and  the  displacements  of 

the  center  of  gravity  of  the  vehicle  from  this  reference  frame  are  VCt) 
and  $(t)  . 

13.  The  velocities  v^^  and  Vy  (relative  to  the  origin  of  the 
coordinate  system)  as  well  as  the  velocities  v^  and  v^  are  related  to 
the  instantaneous  velocity  v  of  the  CG  by 

V  =  |/v^  +  Vy~=  |/vj  +  vj  (10) 

The  side-slip  angle  a  ,  which  is  the  angle  between  the  velocity  vector  v 
and  the  longitudinal  X  axis  of  the  vehicle,  is  related  to  the  velocities 
v„  and  Vy  as 


a 


tan 


_1  ^ 


da 
’  dt 


(11) 


The  yaw  angle  w  and  the  directional  angle  0  are  related  to  a  as 


dO  _  do)  da 

’  dt  dt  dt 


(12) 


Substitution  of  Equation  11  into  Equation  12  leads  to 


11 


dt 


(13) 


14.  The  radius  of  curvature  of  the  trajectory  of  the  center  of  gravity 
(i.e.,  the  distance  between  CR  and  CG,  Figure  3)  is 


3  -  /d0 

1  v/TT  = 

o  /  dt 


2  djo  ^ 

dt  "  '^X  dt  dt 


(14) 


The  coordinates  of  the  trajectory  of  the  center  of  gravity  of  the  vehicle 
can  be  written  as 


'i'(t) 


$(t) 


-P 


V  cos  0  dt 


•'rv 


V  sin  0  dt 


(15) 


15.  The  coordinates  of  the  instantaneous  center  of  rotation  (ICR)  of 
the  hull  in  the  XY  systems  (X^.  ,  Yj.)  and  the  instantaneous  radius  of 
curvature  are  (Figure  3) 


=^1  ■  f  +  S  ■  "y/t  +  <=x  > 

''i-'-'-xf 

Rj  =  n1r2  + 

Track  Slip  Velocity  and  Displacement 

16.  Assume  that  v^j^  (v^j^  =  A^)  and  Vg2  (Vg2  ~  ^2)  are  the  slip 
velocities  of  geometrically  similar  points  of  the  outer  track  and  the  inner 
track,  respectively.  The  X  and  Y  components  of  these  velocities  can  be 
shown  to  be 
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_  _  d(o  _  _  T  d(0 
'^sXl  ■■  ^1  dt  ^1  ^  dt 


"^sYl  “  ^  dt  "  dt  “  "^Y 


For  the  outer  track  (17) 


-  r  ili  -  r  T 
^sX2  ^2  dt  ■  ^2  ^  dt 


For  the  inner  track 


(18) 


’'sY2  ”  "^sYl 


The  angular  velocity  du/dt  and  R  can  be  written  as 


dt  "  bL  ^'^Xl.“  ''sXl  ~  '^X2  ^^5X2^  ^ 


^  ^'^Xl  "  ''sXl  ''X2  "  '^sX2^ 

dt 


(19) 


where  =  the  velocity  of  the  outer  track  relative  to  the  hull 
Vjj2  “  the  velocity  of  the  inner  track  relative  to  the  hull 


The  ratio  of  v^^^  and  v^2  defined  as  the  steering  ratio  e  .  Thus, 


^  “  '^X1^''X2 


Substitution  of  Equations  16  and  20  into  Equation  19  leads  to 


-  (v^  +  For  the  outer  track 


bL  dw 

'^sXl  =  ^'^X2  -  \^X  2 


(20) 


(21) 


^^8X2  =  ''X2  “  ('^X  ■ 


(22) 
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Comparison  between  Equations  21  and  22  and  Equations  17  and  18  results  in 


The  slip  velocities  and  displacements  of  the  outer  and  inner  tracks  can  be 
obtained  from  Equations  17,  18,  21,  and  22,  Thus, 


where  t^^  =  (L/2  - 

t2  =  (L/2  -  X)/Vjj2 

=  initial  displacement  of  the  outer  track 
Aj2  =  initial  displacement  of  the  inner  track 

The  balance  of  forces  and  moments  dictates  that  these  initial  displacements 
be  numerically  equal  to  Lj{  (jj  is  the  coefficient  of  rolling  resistance 
which  must  be  measured  experimentally  or  calculated  from  empirical  relations 
such  as  those  given  in  Reference  8  for  each  soil  type  and  each  vehicle). 
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Inertial  Forces 


17.  The  X  and  Y  components  of  the  inertial  force  can  be  shown  to  be 
(Reference  7) 


F 


CX 


w 

g 


+  V 


doA 

Y  dt/’ 


CY 


W  /^^Y  doA 
g  \dt  "  '^X  dt/ 


The  Rolling  Resistance 


(28) 


18.  The  rolling  resistance  is  a  function  of  terrain  type,  vehicle 
speed,  track  condition,  etc.  Therefore,  rolling  resistance  should  be 
measured  for  every  specific  condition.  In  this  formulation,  however,  the 
rolling  resistance  is  assumed  to  be  proportional  to  normal  load.  Thus, 

1 

Rg  =  i  j  ^  [rj^(x)  +  r2(x)]dx  (29) 

2 


Equations  of  Motion 

19.  Steerability  and  stability  of  tracked  vehicles  depend  on  the 
dynamic  balance  between  all  forces  and  moments  applied  on  the  vehicle. 
According  to  Figure  4,  the  following  three  equations  govern  the  motion  of 
the  vehicle; 


Etj^(x)  + 


t2(x)]  dx  - 


[rj^(x)  +  r2(x)]  dx  =  f^^ 

a 


+  92  (x)]  dx  =  fj^Y 


(30) 


(31) 
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1  1 
j  ^  [q^^Cx)  +  q2(x)]  (x  -  c^)  dx  +  jJ  [tj^(x) 


t„(x)]dx 


1 

’2 


I_  ,2 


dt 


(32) 


where  tj^(x)  =  dL^  Tj^(x)/W  ,  t2(x)  =  dL^  :j:2(x)/W  ,  q^^Cx)  =  dL^  Qj^(x)/W  , 
q2(x)  =  dL^  Q2(x)/W  ,  =  F^^/W  ,  and  =  F^^/W 


and  I  =  mass  moment  of  inertia  about  an  axis  passing  through  the  center  of 
2 

gravity  of  the  vehicle  and  parallel  to  the  Z  axis  (Figure  3) .  Equations 
30  through  32  with  the  aid  of  Equations  7  through  29  constitute  three  equa¬ 
tions  that  Involve  three  unknowns.  The  three  unknowns  are  either  Vj^  , 
v^  ,  and  dw/dt  or  ,  §2  >  P  •  order  to  obtain  a  complete  solu¬ 

tion  for  either  of  the  two  sets  of  unknowns,  one  of  the  following  driving 
conditions  must  be  specified:  (a)  time  history  of  the  steering  ratio  e(t) 
and  the  initial  speed  of  the  vehicle,  (b)  time  history  of  the  velocity  of 
the  individual  tracks  and  Vj^(t)  and  the  initial  speed  of  the 

vehicle,  (c)  time  history  of  the  velocity  of  the  vehicle  v(t)  and  the 
trajectory  of  motion,  (d)  time  history  of  the  velocity  of  the  vehicle  and  a 
constant  value  of  steering  ratio  e  ,  or  (e)  the  trajectory  of  motion  and  a 
determination  of  the  maximum  velocity-time  history  at  which  the  vehicle  can 
traverse  the  specified  trajectory.  A  computer  program  called  AGIL  was 
developed  to  solve  Equations  30  through  32  using  Newton's  Iteration  tech¬ 
nique.  In  addition,  this  computer  program  has  the  capability  of  calculating 
the  total  power  requirements  (PT)  as  well  as  the  power  required  at  each 
sprocket  (Reference  7). 


PART  IV:  CORRELATION  OF  THEORETICAL  PREDICTIONS.  WITH  TEST  DATA 


20.  In  order  to  determine  the  accuracy  and  range  of  application  of 
the  terrain-vehicle  model  a  series  of  slalom  and  circular  turn  tests  was 
conducted,  in  several  areas  in  the  vicinity  of  Vicksburg,  Mississippi..  The 
results  of  these  tests  are  presently  being  analyzed  for  correlation  and’ 
comparison  with  the  model  predictions.  Data  from  five  turn  tests  (test, 
series  107  through  111),  however,  have  been  reduced  and  are  available  for 
correlation  analysis.  These  data  offer  the  first  opportunity  to  evaluate 
the  accuracy  of  the  terrain-vehicle  model.  A  brief  description  of  the  test 
conditions  for  this  series  of  tests  and  the  analysis  of  the  data  and  com¬ 
parison  with  the  model  predictions  are  presented  in  the  subsequent  sections. 

Experimental  Program 

21.  Test  series  107  through  111  was  conducted  in  an  area  termed  as 
New  Ground  Area,  located  north  of  Redwood,  Mississippi.  The  area  is  barren 
and  iinplowed.  The  soil  at  this  area  is  a  soft  plastic  clay  classified’  as  CH 
according  to  the  Unified  Soil  Classification  System.*  Each  test  involved 
the  steering  of  the  vehicle  in  a  circular  path  by  first  accelerating  the 
vehicle  to  maximum  speed  (controlled  by  either  the  available  power  or  the 
actual  physical  stability  of  the  vehicle)  and  then  continue  turning  with  a 
more  or  less  constant  speed.  The  tracked  vehicle  used  for  these  experiments 
is  an  armored  personnel  vehicle  with  the  characteristics  given  in  Table  1, 
The  actual  data  collected  during  each  test  consisted  of  time  histories  of 
(a)  the  inner  and  outer  track  velocities,  (b)  the  speed  of  the  vehicle,  (c) 
the  turning  radius,  and  (d)  the  power  requirement.  The  lateral  acceleration 
of  the  vehicle  was  then  calculated  from  items  b  and  c  above.  In  addition  to 
the  above  data,  numerous  slow  and  fast  in  situ  direct  shear  tests  were 
conducted  at  each  test  site  in  order  to  characterize  the  soil  in  terms  of 
the  soil  parameters  given  in  Equations  1  through  4.  Table  2  summarizes  the 
soil  parameters  for  the  New  Ground  Area.  Also  shown  in  Table  2  is  the 
coefficient  of  rolling  resistance  for  the  vehicle  which  was  obtained  experi¬ 
mentally  (through  torque  measurements  for  straightforward  motion)  at  the 
site. 

*  The  Unified  Soil  Classification  System  is  described  in  Reference  9. 
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Table  2 

Suiamary  of  the  Soil  Parameters  for  the  New  Ground  Area 

"n5tial  shear  stiffness  coefficient  (G)  =  200  psi/in. 

Cohesive  strength  of  the  material  corresponding  to 

static  loading  (C  *  A  -  M)  =0.94  psi 

Static  failure  envelope  (Equation  1)  parameter  (A)  =  5.00  psi 

Static  failure  envelope  (Equation  1)  parameter  (M)  »  4.06  psi 

Static  failure  envelope  (Equation  1)  parameter  (N)  ■  0.22  1/psi 

Dynamic  fa^ure  envelope  (Equation  2)  parameter  (C^)  «  0.61  psi 

Dynamic  failure  envelope  (Equation  2)  parameter  (A)  «  3.68  sec/ln. 

Coefficient  of  rolling  resistance  (j()  ■  0.2 

WES  cone  index  (Cl)  for  0-  to  6-ln.  depth  «  30  psi 

WES  cone  index  (Cl)  for  6-  to  12-in.  depth  ■  65  psi 


The  soil  parameters  A,  M,  and  N  in  Table  2  were  determined  by  fitting 
Equation  1  to  the  experimental  static  failure  envelope  of  the  soil  (from 
slow  tests) .  The  initial  shear  stiffness  modulus  G  was  then  determined’  by 
fitting  Equations  1  and  3  to  the  experimental  shear  stress-shear  deformation 
curves  from  slow  tests.  Finally  the  parameters  and  A  were  obtained  by 
fitting  Equations  2  and  4  to  the  experimental  failure  envelope  and  shear 
stress-shear  deformation  curves  from  fast  tests.  The  representative  failure 
envelope  obtained  from  experimental  static  data  and  the  corresponding 
analytical  curve  (Equation  1)  are  shown  in  ‘Figure  5.  It  is  observed  from 
Figure  5  that  the  agreement  between  the  two  curves  is  excellent.  Typical 
experimental  and  the  corresponding  analytical  static  and  dynamic  shearing 
stress-shearing  deformation  curves  for  the  clay  soil  from  the  New  Ground 
Area  are  shown  in  Figure  6.  Again,  the  experimental  results  and  the  model 
simulations  are  in  close  agreement. 

Theoretical  Predictions 


22.  The  purpose  of  the  circular  turn  tests  was  to  create  a  condition 
of  steady-state  motion  for  the  vehicle  to  evaluate  the  accuracy  of  the  model 
for  this  simple  mode  of  motion.  Unfortunately,  due  to  factors  such  as 
surface  roughness,  inhomogeneity  of  surface  materials,  and  the  fact  that  the 
ground  is  not  an  ideally  flat,  level  surface,  it  is  not  possible  to  maintain 
the  vehicle  in  a  perfect  steady-state  mode  of  motion.  These  factors  adverse¬ 
ly  affect  the  ability  of  the  driver  to  steer  the  vehicle  at  a  constant 
radius  with  a  constant  velocity.  Therefore,  the  motion  of  the  vehicle 
during  the  entire  test  event  cannot  be  considered  as  steady-state  motion. 

For  this  reason,  for  each  test  a  time  window  was  selected  where  the  motion 
of  the  vehicle  could  be  reasonably  approximated  as  steady  state.  The  time 

windows  and  the  corresponding  test  data  consisting  of  the  steering  ratio  e  , 

2 

vehicle  velocity  v  ,  turning  radius  ,  lateral  acceleration  v  /R^B 
total  power  FT  are  given  in  Table  3  for  all  the  tests.  The  computer 
program  AGIL  was  used  to  simulate  these  steady-state  test  conditions  using 
the  vehicle  characteristics  and  the  soil  parameters  from  Tables  1  and  2, 
respectively.  Comparisons  of  the  model  predictions  with  experimental  data 
are  presented  in  Figures  7  through  12  for  turning  radius  versus  steering 
ratio,  inner  track  velocity  versus  turning  radius,  outer  track  velocity 


21 


/ 


Tab  le  3 


Steadprbuate.  Time  VEndows  and  Corresponding  Test  Date 


Teat 

Nbv 

sec 

Steerlii$t 

Ratio.' 

e. 

TeMcle 
Speedi 
V.  mpbi 

107 

60)  6& 

1.08! 

24.88! 

108! 

00 

% 

1.29 

14.9-31 

109) 

66 

1.12 

1. 10! 

15;.73. 

18'.22= 

iia 

55.  5.9' 

1.26 

16.08'. 

111 

48'  51 

1.57 

11.79 

Turning, 

Radius 

LateraL 

Acceleration; 

v:  /•  V' 

Fow^ 
•Fir,  Iffi' 

155.77 

O'.  265 

210 

6.3.69! 

a.234 

159 

103'.  86 

0).16. 

177 

121.26 

0.18 

184 

83.40- 

0..21 

200) 

33.74 

0.275. 

211 

22 


- FIELD  MEASUREMENTS 

-  MODEL  BEHAVIOR 


SHEAR  STRESS  T,  PSI  SHEAR  STRESS  T  PSI 


LEGEND 


SHEAR  DISPLACEMENT  A  ,  INCH 


SHEAR  DISPLACEMENT  A  ,  INCH 


Figure  6.  Comparison  of  experimental  shearing  stress-shearing 
displacement  relation  with  model  behavior 
(o  -•  4.14  psi)  for  clay  soil  from  New  Ground  Area 
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model  predictions  with  experimental  data 


FIELD  MEASUREMENTS 
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predictions  with  experimental  data 
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model  predictions  with  experimental  data 


versus  turning  radius,  vehicle  speed  versus  turning  radius,  power  require¬ 
ment  versus  turning  radius,  and  lateral  acceleration  versus  turning  radius, 
respectively.  The  model  predictions  in  Figures  8  through  12  are  based  on 
both  power  cutoff  and  preliminary  stability  criteria.  The  turning  radius¬ 
steering  ratio  relation  showri  in  Figure  7,  however,  is  unique  for  a  given 
vehicle  and  soil  condition.  The  power  cutoff,  as  indicated,  is  controlled 
by  the  available  power.  According  to  Reference  7,  the  preliminary  stability 
criteria  are  based  on; 

a.  Rapid  change  in  the  slip  velocity  of  the  inner  or  the  outer 
track. 

_b.  The  pivot  point  moves  outside  the  front  edge  of  the  track- 
ground  contact  area  (i.e, ,  the  offset  equals.  0.5  L  when  the 
center  of  gravity  and  center  of  geometry  of  the  vehicle 
coincide) . 

£.  Rapid  decrease  or  increase  in  the  turning  radius. 

These  stability  conditions  usually  take  place  at  different  vehicle  veloci¬ 
ties,  The  unstable  vehicle  velocity  is  chosen  as  the  minimum  of  these 
velocities.  For  comparison  with  the  experimental  data  the  lower  vehicle 
velocity  corresponding  to  either  the  stability  criteria  or  the  power  cutoff 
condition  must  be  selected.  As  Indicated  in  Figure  10,  for  the  turning 
radii  of  34  ft  (test  111)  and  64  ft  (test  108)  stability  criteria  control 
the  velocity  of  the  vehicle.  For  the  turning  radii  of  83  ft  (test  110), 

104  ft  and  121  ft  (test  109),  and  156  ft  (test  107)  the  velocity  of  the 
vehicle  is  controlled  by  the  available  power.  With  this  in  mind,  the 
experimental  data  in  Figures  7  through  12  compare  very  favorably  with  the 
corresponding  model  predictions.  This  is  particularly  true  in  the  case  of 
track  velocities  and  vehicle  speed  (Figures  8  through  10).  Slight  observable 
differences  between  the  data  and  model  predictions  in  Figures  7  through  12 
should  be  expected  because  of  the  deviations  in  the  test  conditions  from  the 
steady-state  mode  of  motion. 

23,  To  demonstrate  the  ability  of  the  terrain-vehicle  model  for 
treating  transient  motion,  the  response  of  the  vehicle  for  test  107  was 
simulated  for  the  entire  test  event.  The  measured  time  histories  of  the 
inner  and  outer  track  velocities  were  used  to  drive  the  model.  For  these 
specified  driving  conditions,  the  time  histories  of  the  vehicle  speed,  power 
requirements  and  lateral  acceleration,  and  the  trajectory  of  the  center  of 
gravity  of  the  vehicle  were  then  predicted  and  compared  with  the  correspond¬ 
ing  field  measurements.  Figure  13  depicts  the  time  histories  of  the  inner 
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Figure  13.  Outer  and  inner  track  velocity-time  histories  for 
test  107;  field  measurement  and  filtered  data 
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and  outer  track  velocities.  The  actual  field  measurements  are  quite  noisy, 
particularly  for  times  greater  than  approximately  60  sec.  These  high- 
frequency  oscillations  are  believed  to  be  mostly  due  to  instrumentation  and 
were  filtered'  out.  The  filtered  records  are  also  shown  in  Figure  13  and  are 
simply  "best"  fit  curves  to  the  field< -measurements  satisfying  the  condition 
that  the  total  area  under  both  curves  should  be  equal.  These  filtered  track 
velocity-time  histories  were  used  as  input  to  drive  the  terrain-vehicle 
model.  Comparisons  of  the  predicted  time  histories  of  the  vehicle  speed, 
power  requirements,  and  lateral  acceleration  with  the  corresponding  field 
measurements  are  shown  in  Figures  14  through  16,  respectively.  Similar  to 
Figure  13,  the  field  measurements  are  quite  noisy.  As  anticipated,  the 
predicted  results  do  not  manifest  these  oscillations  because  of  the  filter¬ 
ing  of  the  input  data.  The  degree  of  correlation  of  the  predicted  and 
measured  results,  however,  is  quite  good,  indicating  that  the  modeling  of 
the  overall  interaction  between  the  terrain  and  the  track  is  physically 
reasonable.  Comparison  of  the  predicted  and  measured  trajectory  of  the 
center  of  gravity  of  the  vehicle  is  shown  in  Figure  17.  Again,  the  com¬ 
parison  between  the  predicted  and  measured  results  is  quite  favorable. 
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Figure  15.  Total  power-time  history  for  test  107 
comparison  of  model  predictions  with 
experimental  data 


PART  V:  SUIffi^Y  AND  CONCLUSION 

24.  A  mathematical  model  of  ter rain- vehicle  interaction  for  predict¬ 
ing  the  steering  performance  of  track-laying  vehicles  has  been  developed  and 
computerized  for  numerical  application.  The  model  contains  some  of  the 
basic  parameters  governing  the  steering  performance  of  tanks,  such  as  track 
slippage,  centrifugal  forces,  vehicle  characteristics,  and  soil  type.  The 
model  is  applicable  to  both  steady-state  motion  and.  trahsient  motion  under 
changing  control  commands.  The  model  is  partially  validated  by  comparing 
field  measurements  for  a  specific  vehicle  with  the  corresponding  model 
predictions.  Efforts  are  presently  under  way  at  WES  to  compare  the  model 
predictions  with  experimental  data  for  a  broader  range  of  test  conditions. 
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APPENDIX 'A 
NOTATION 

a  AL^/W 

A  Soil  parameter 

b  B/L  ; 

B  Track  tread 

c  CL^/W 

c'  C,L^/W 

d  d 

C^/L 

C  Static  cohesive  component  of  shear  strength 

Added  cohesive  strength  due  to  dynamic  loading 

Cjj  Abscissa  of  the  center  of  gravity  of  the  vehicle 

Slip  radius  of  the  outer  track 

C2  Slip  radius  of  the  inner  track 

CG  Center  of  gravity  of  the  vehicle 

Cl  WES  cone  index 

CR  Center  of  rotation  of  the  vehicle 

d  D/L 

D  Track  width 

^cx 

*CY 

Fg  Inertial  force 

F  ^  Longitudinal  component  of  inertial  force 

F^Y  Transverse  component  of  inertial  force 

^  Coefficient  of  rolling  resistance 

g  Acceleration  due  to  gravity 

G  Initial  shear  stiffness  coefficient 

h  H/L 

K  Height  of  center  of  gravity 

I  Mass  moment  of  inertia  of  the  vehicle  about  an  axis  passing  through 

its  center  of  gravity  and  parallel  to  the  Z  axis 

ICj^  Center  of  slip  rotation  of  the  outer  track 
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IC2  Center  of  slip  rotation  of  the  inner  trac,c 

ICR  Instantaneous  center  of  rotation  of  the  vehicle 
L  Contact  length  of  track 

m  ML^/W 

M  Soil  parameter 

n  NW/L^ 

N  Soil  parameter 

0  Center  of  geometry  of  the  vehicle 
p  P/L 

P  Offset  (distance  from  center  of  gravity  to  pivot  point  of  vehicle) 
PT  Total  power  =  PTl  +  PT2 

PTl  Power  required  by  the  sprocket  of  the  outer  track 

PT2  Power  required  by  the  sprocket  of  the  inner  track 

PTD  Differential  power  =  PTl  -  PT2 

qj^(x)  dL^q^(x)/W 

q2(x)  dL^Q2(x)/W 

Q^(X)  Transverse  component  of  shear  stress  along  the  outer  track 

Q2(X)  Transverse  component  of  shear  stress  along  the  inner  track 

rj^(x)  dL^Rj^(x)/W 

r2 (x)  dL^R2 (x) /W 

R  Ordinate  of  the  instantaneous  center  of  rotation  of  the  vehicle 
R^  Radius  of  the  trajectory  of  the  center  of  gravity  of  the  /ehicle 

Rg  Rolling  resistance 

Rj^(X)  Normal  stress  under  the  outer  track 

R2(X)  Normal  stress  under  the  inner  track 

Rj.  Instantaneous  radius  of  curvature 

t  Time 

t^(x)  dL^Tj^(x)/W 

t2(x)  dL^T2(x)/W 
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Tj^(X)  Longitudinal  component  of  shear  stress  along  the  outer  track 

T2(X)  Longitudinal  component  of  shear  stress  along  the  inner  track 

V  Velocity  of  the  vehicle 

Total  slip  velocity  of  the  outer  track 

Vg2  Total  slip  velocity  of  the  inner  track 

Vsxi  Longitudinal  component  of  slip  velocity  of  the  outer  track 

VsX2  Longitudinal  component  of  slip  velocity  of  the  inner  track 

^sYl  Transverse  component  of  slip  velocity  of  the  outer  track 

v^Y2  Transverse  component  of  slip  velocity  of  the  inner  track 

v^  Longitudinal  component  of  velocity  of  the  vehicle 

v^j^  Longitudinal  component  of  velocity  of  the  outer  track 

v^2  Longitudinal  component  of  velocity  of  the  inner  track 

v^  Transverse  component  of  velocity  of  the  vehicle 

W  Weight  of  the  vehicle 

X  X/L 

X,Y,Z  Local  coordinate  system 
y  Y/L 

2  Z/L 

a  Side-slip  angle 

Angle  of  slip  direction  of  the  outer  track 

Y2  Angle  of  slip  direction  of  the  inner  track 

A  Shearing  deformation 

Ajj^  Initial  displacement  of  the  outer  track 

Aj.2  Initial  displacement  of  the  inner  track 

Aj^  Shearing  deformation  of  soil  under  the  outer  track 

Aj^  Time  rate  of  shearing  deformation 

A2  Shearing  deformation  of  soil  under  the  inner  track 

^2  Time  rate  of  shearing  deformation 
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h  *2''- 

e  Steering  ratio 

0  Directional  angle 

A  Material  constant  related  to  rate  effect 

A  AL 

y  GL^/W 

5l  C^/L 

?2  C^/L 

a  Normal  stress 

T  Shear  stress 

Tj^  Maximum  shear  strength 

Coordinate  system  fixed  on  level  ground 
u)  Yaw  angle 


A4 


In  accordance  with  letter  fron  bAEX'-RDG,,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsiraile  Catalog  Cards,  for 
Laboratory  Te.chnical  Publications,  a  facsimile;  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Baladi,  George  Y. 

Analysis  of  steerability  of  tracked  vehicles; 
theoretical  predictions  -versus-  f ieldTroeasuremenijs  : 

Final  report  /  by  George  Y.  Baladi,  Behza”d "Sohani 
(Structures  Laboratory,.  U.S.  Army  Engineer  Waterways 
Experiment  Station)  ;  prepared  for  Office,  Chief  of 
Engineers,  U.S.  Army  ;  monitored  by  Geotechnical 
Laboratory,  U.S.  Army  Engineer  Waterways  Experiment 
Station  ;  Springfield,  Va.  :  available  from  NTIS, 

1981. 

39,  1*  ;o.  :  ill.  ;  27  cm.  —  (Miscellaneous  paper  /  U.S. 
Army  Engineer  Waterways  Experiment  Station  ;  SL-81-3) 
Cover  title. 

"March  I98I." 

"Under  Project  liAl6ll02AT22  Task  CO,  Work  Unit  001." 
Bibliography:  p.  39. 

1.  AGIL  (Computer  program).  2.  Computer  programs, 

3.  Mathematical  models.  L,  Soil  mechanics.  5.  Track¬ 
laying  vehicles.  6.  Vehicles,  Military.  I.  Bphani, 
Behzad.  II.  United  States,  Army.  Corps  of  Engineers. 


Baladi,  Gec  rge  I. 

Analysis  of  steerability  of  tracked  vehicles  :  ...  1981. 

Office  of  the  Chief  of  Engineers.  III.  United  States. 
Array  Engineer  Waterways  Experiment  Station.  Structures 
Laboratory.  IV.  Title  V.  Series:  Miscellaneous 
paper  (United  States.  Army  Engineer  Waterways  Experiment 
Station)  ;  SL-81-3. 

TA7.W3lim  no.  SL-81-3 


